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Abstract 

We simulate the transformation of a classical fluid into a quantum-like (super)- 
fluid by the application of a generalized quantum potential through a retro- 
active loop. This numerical experiment is exemplified in the case of a non- 
spreading oscillating wave packet in a harmonic potential. We find signatures 
of a quantum-like behavior which are stable against various perturbations. 

1 Introduction 

One of us has recently proposed a new general concept of macroscopic quantum-type 
laboratory experiments [1]. It consists of applying, through a real time retroactive loop, a 
generalized quantum potential on a classical system. Indeed, one can show that the system 
of equations (Euler equation and continuity equation) that describes a fluid in irrotational 
motion subjected to such a generalized quantum potential, that reads Q = —2V 2 A y /p/y r p 
in terms of the density p, is equivalent to a generalized Schrodinger equation. In this 
derivation, the quantum potential is no longer founded on the quantum Planck's constant 
h, but on a new constant V which can take any macroscopic value. While it would be 
impossible with present days technology to simulate standard quantum effects by this 
method because of the smallness of P, which is given by V = h/2m in standard quantum 
mechanics, the use of a macroscopic value for this constant nevertheless preserves some 
of the properties of a quantum-like system. Namely, its density distribution is given 
by the square of the modulus of a complex function which is solution of a Schrodinger 
equation. Therefore such a system is expected to exhibit some quantum-type, superfluid- 
like macroscopic properties (though certainly not every aspects of a genuine quantum 
system) . 

In the present paper we validate this concept by numerical simulations of a fluid 
subjected to such a generalized quantum force, as an anticipation of a future real labora- 
tory experiment. The example chosen for this first attempt is the appearance of a non- 
spreading quantum-like oscillating wave packet in a compressible fluid (e.g., a plasma) 
subjected to an attractive harmonic oscillator potential. 



2 Theoretical background 



We consider a classical macroscopic compressible fluid described by the Euler and the 
continuity equations: 

l + v-vjv = -V<f>, (1) 

^ + dw(pV) = 0, (2) 

where (p is an exterior scalar potential. We assume as a first step that the pressure term 
is negligible and that the fluid motion is potential, i.e., 

V = VS. (3) 

We now assume that we apply to the fluid (using density measurements and a retroaction 
loop) a varying force which is a function of the fluid density in real time, namely, a 
"quantum-like" force Fq deriving from the potential 

Q = _ 2V 2 / Wl (4) 

This potential is a generalisation of the standard quantum potential [2], since here the 
constant V can have any value, while in standard quantum mechanics it is restricted to 
the only value V = %/2m. As recalled in what follows, this generalization still allows to 
recover a Schrodinger-like equation. 

The Euler and continuity system becomes 

| + V.V)V = -V(«-^), (5, 

g+div(^) = 0, (6) 

The system of equations (5,6) can then be integrated under the form of a generalized 
Schrodinger equation. 

Indeed, equation (5) takes the successive forms 

^(VS) + iv(VS) 2 + vL-2D 2 ^j =0, (7) 

V (f + ±(VS)^- 2I >^)=0, (8, 
which can be integrated as 

f + l(VSf + (j) + K-2V^ = 0, (9) 

where K is a constant that can be renormalized by a redefinition of the potential energy 
(f). Let us now combine this equation with the continuity equation as follows: 



Finally we set 

rf, = VP x e^, (11) 
and the equation (10) is strictly identical to the following generalized Schrodinger equation 

D 2 AV + iV^ - ^ = 0, (12) 

as can be checked by replacing in it ip by its expression (11). Recall that such an equation 
has also been directly obtained, in terms of a density of probability instead of a density 
of matter, as the integral of the equations of geodesies in a nondifferentiable space-time 
[3, 4]. Given the linearity of the equation obtained, one can normalize the modulus of 
tp by replacing the matter density p by a probability density P — p/M , where M is the 
total mass of the fluid in the volume considered: this will be equivalent. 

The solutions ip = \ip\ x exp(z#) of this equation directly provide the density and the 
velocity field of the fluid at every point, namely 

V = 2VV6, p = M\^\ 2 . (13) 

Its imaginary part and its real part amount, respectively, to the continuity equation, and 
to the energy equation that writes: 

dS 1 o 9 Av/p , 

The above transformation from the fluid mechanics-like equations to the Schrodinger- 
type equation is similar to a Madelung transformation [5], but it is here performed in the 
reversed way and generalized to a constant different from ft /2m. 

It could be therefore possible by this method to simulate a "Schrodingerian" system, 
e.g., a partly quantum-like superfluid system coming under two of the axioms of quantum 
mechanics, namely, it is described by a wave function ip such that p oc \ip\ 2 , which is 
solution of a Schrodinger-type equation. 



3 Application to the oscillating wave packet 

As an example of application and as a preparation for a laboratory experiment, let us con- 
sider the simplified case of one-dimensional fluid motion in an external harmonic oscillator 
potential (p = (l/2)u 2 x 2 . This system is described by the two following equations: 

Here we have written the continuity equation in terms of lnp, which will be useful in 
the numerical simulations that follow. These two equations are equivalent to the one- 
dimensional generalized Schrodinger equation: 



It is well known that it is possible to find a solution of this equation in the form of a wave 
packet whose center of gravity oscillates with the period of the classical motion and which 
shows no spreading with time [6, 7, 8]. Assuming that the maximal amount by which the 
center of gravity is displaced is a, the wave function reads in this case 



i 



-0 = ^_^_^ 4 e ~^(x-a cos u;t) 2 x e -i(\ujt+^ax sin^-^a 2 sin 2ut) _ 

Therefore the probability density reads 



p = |^|2 = e -^(x-accsujt)\ ( 19 ) 

This is an interesting case for a test of a genuine quantum behavior, since it involves a 
non vanishing phase in an essential way although this is a one-dimensional system. The 
velocity field is given by 

V = —ausinut, (20) 
while the expression for the quantum potential is 

Q(x, t) = Vcj — ^j?{x — acosut) 2 , (21) 
so that the quantum force writes 

dO 

F Q cj 2 (x - acosujt). (22) 

Therefore the (varying) energy takes the form 

E = -V 2 + (f) + Q = Vcj + au 2 x cos out - -a 2 u 2 cos(2cjt). (23) 

When it is applied to the center of the packet x = acosut, this expression becomes 

E c = Vcj + ^a 2 u; 2 . (24) 

We recognize in the second term, as expected, the energy of a classical pendulum. Con- 
cerning the first term, since standard quantum mechanics corresponds to the particular 
choice V = h/2m (here with m = 1), the term Vuj is the generalization of the vacuum 
energy for an harmonic oscillator, ^Ttuo. 

Therefore we verify that the application of a quantum potential on the fluid has given 
to it some new properties of a quantum-like nature, such as a zero-point energy and the 
conservation of the shape of the wave packet. 



4 Proposed laboratory experiment 

In order to prepare a real laboratory experiment aiming at achieving such a new macro- 
scopic quantum-like (super)fluid, we shall now present the result of numerical simulations 
of such an experiment. To this purpose these simulations are not based on the Schrodinger 
form of the equations, but instead on the classical Euler + continuity equations and on 
the application by feedback of a generalized quantum-like force. 



The suggested experiment consists of: 

(i) measuring with detectors the density at regular time interval {t n } on a grid at 
positions {xj}; 

(ii) computing from these measurements the quantum force (Fq) u = 2V 2 V(Ay / p^/ y/fa) 
at each time t n \ 

(iii) applying the new value of the force to the fluid at each time t n , therefore simulating 
by such a feedback the presence of a quantum-like potential. 

The advantage of such a proposal is that one is no longer constrained by the standard 
quantum value V = h/2m that fixes the amplitude of the quantum force, and that one 
can therefore give to it a macroscopic value, vary it, study its transition to zero (quantum 
to classical transition), etc... 



5 Iterative fitting simulation 

In this first simulation, we assume that the quantum force (which is a third derivative of 
the density) is not computed directly from the values of the density, but from a polynomial 
fit of the distribution of In p. In the special case considered here (the oscillating wave 
packet), we use a Gaussian fit of the density distribution (i.e. a second order polynomial 
fit to In p), so that we need to know only the mean and dispersion. More generally, one can 
decompose the distribution of In p(x) into its successive moments. Therefore the density 
is written as 



p n (x) oc exp 



1 / X x n x 



(25) 



so that, once the mean and dispersion x n and a n at time t n are computed, the quantum 
force to be applied at each step (n) writes: 

(FqUx) = p2 °V n) - (26) 



5.1 Numerical simulation 

Our numerical simulation is performed by a simple Mathematica program which repro- 
duces the steps of the real experiment, namely, at each time step t n : 

(i) We compute the mean and the dispersion of positions x according to the density 
distribution: 

x = ^pix^Xj/^pixj), (27) 

3 3 

a 2 = '£ P (x j )(x j ~ x) 2 /^2p(xj). (28) 

3 3 

(ii) The force Fq to be added then writes in terms of these quantities 
where Sx is the grid interval and intervenes here because we use finite differences. 



(iii) We compute the logarithm of the density In p and the velocity V at next time step 
t n+ l by transforming equations (15, 16) into centered finite-difference equations (Forward 
Time Centered Space, FTCS scheme) using the Lax-Friedrichs method [9], namely, 

. , In p 7 ! , -, + In p 7 } t St ( -\ 

yn I yn / _ T/n \ 

v?+ 1 = ^4 - + st {- v ? J+1 2fa J " 1 + + • (3D 

The lower index (j) is for space x and the upper one (n) is for time t; 5t is the time step 
and F(x) — —uo 2 x is the external harmonic oscillator force. In the above Lax method, 
the terms In p™ and VJ 1 are replaced by their space average, which has the advantage to 
stabilize the FTCS scheme. 

The initial conditions are given by the density distribution (Eq. 19) for t = 0. 

Although this is a simple scheme (we have not attempted at this stage to better control 
numerical error diffusion), it has given very encouraging results, since it has reproduced on 
several periods the expected motion of the quantum oscillating wave packet (see Figure 1). 




Figure 1 : Result of the numerical integration of a Euler + continuity one dimensional system with generalized 
quantum potential for the oscillating wave packet in an harmonic oscillator field. The quantum force applied on 
the fluid is calculated from a gaussian fit of the density distribution. The figure gives the density distribution 
obtained in function of position (space grid from 25 to 75) and time (time steps from 1 to 75, i.e., 1.2 period). 



5.2 Perturbation of initial conditions 

One of the possible shortcomings in the passage from the simulation to a real experiment 
may come from fluctuations in the initial conditions. Indeed, in the previous simulations, 
we have taken as initial density distribution that of the exact quantum wave packet. In 
order to be closer to a real experimental situation, we have therefore performed a new 
simulation similar to that of Sec. 5.1, but with an initial density distribution that is 
perturbed with respect to the Gaussian solution (Eq. 19): we have multiplied its values 
p(xj) at each point {xj} of the space grid by exp(a), where a is random in the interval 



[0,1]. A typical resulting initial density distribution is given in Figure 2, followed by 
the distributions obtained on a full period (sub-figures 1 to 12) after application of the 
generalized quantum force. 

Once again the result obtained is very enrouraging as concerns the possibility of per- 
forming a real laboratory experiment, since, despite the initial deformation, the wave 
packet remains stable during several periods. Moreover, not only the mean and dispersion 
of the evolving density distribution remain close to the ones expected for the quantum 
wave packet, but, as can be seen in Figure 2, the initial perturbations have even been 
smoothed out during the feedback process. 



5.3 General account of uncertainties 

This encouraging result leads us to attempt a numerical simulation under far more difficult 
conditions: in order to simulate the various uncertainties and errors that may occur in a 
real experiment, in particular as concerns the density measurement, the application of the 
force, and physical effects not accounted in the simulation such as pressure (see below), 
vorticity, etc., we have now added a fluctuation at each step of the retroactive loop. 
Namely, at each time step t n , we have multiplied the density p(xj) at each point {xj} of 
the space grid by exp(a), where a is random in the interval [0,1]. 

As can be seen in Figure 3, despite the large errors added, the numerical simulation 
shows an oscillating wave packet which, despite its large fluctuations, keeps its coherence. 
In particular, it keeps the values of the mean and dispersion (to about 5 percent) expected 
for the quantum solution on about 1/3 of period before the end of the simulation due to 
numerical errors. 



5.4 Account of pressure 

The addition of a pressure term in the initial Euler equation still allows one to obtain a 
Schrodinger-like equation in the general case when Vp/p is a gradient, i.e., Vp/p = Vw. 
This is the case of an isentropic fluid, and, more generally, of every cases when there is an 
univocal link between pressure and density, e.g., a state equation [10]. The Euler equation 
with quantum potential and external potential reads 

(| + V.v)v=-v(. + TO -2^), (32) 

and it can therefore, in combination with the continuity equation, be integrated in terms 
of a Schrodinger-like equation, 

V 2 A^j + iV^j) - V = 0. (33) 

Now the pressure term needs to be specified through a state equation, which can be 
chosen as taking the general form p — kpp 1 '. The special case 7=1 can be recovered 
and its amplitude established by taking the acoustic limit p = po + p' p — po + p f and 
p' — (? s p j where c s is the sound velocity in the fluid. Therefore one obtains a linear relation 
p = a + clp, so that the pressure term in the Euler equation finally reads Vp/p = fc^Vlnp, 
while w = kplnp = A;^ In | ^ | 2 , with k p = c 2 s . This means that the integrated equation is 
now a nonlinear Schrodinger equation, 

d 1 
V 2 A^ + iV— ip-kphi\ip\ ip = -<f>ip. (34) 
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Figure 2: Result of the numerical integration of a Euler + continuity one dimensional system with added 
quantum potential for the oscillating wave packet. The conditions are the same as in Figure 1, except for the 
addition of a perturbation on the initial density distribution (left top figure). The quantum force applied on 
the fluid is calculated from a Gaussian fit of the density distribution. The successive figures give the density 
distribution obtained in function of position (space grid from 25 to 75) and time (64 time steps corresponding 
to one period, among which twelve of them, equally distributed, are shown). 



In the highly compressible case the dominant pressure term is rather oc p 2 , and the In \ip\ 
term is replaced by \ip\ 2 in the non-linear Schrodinger equation (see e.g. [11]). 
The numerical integration is now performed by generalizing Eq. (31) as 

V i ~ 2 + St {- V i 25~x + F J+( F Qh- k p ^ • 

(35) 

The result is given in Figures 4 and 5 for two different values of the pressure amplitude k p . 



Fi gure 3: Result of the numerical integration of a Euler + continuity one dimensional system with added 
generalized quantum potential for the oscillating wave packet. The conditions are the same as in Figure 1, 
except for the addition of a perturbation on the density distribution at each time step of the simulation. The 
quantum force applied on the fluid is calculated from a Gaussian fit of the density distribution. The successive 
figures give the density distribution obtained in function of position (space grid from 65 to 105) and time (time 
steps 1 to 25, which corresponds to almost half a period). 
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Figure 4: Result of the numerical integration of a Euler and continuity one-dimensional system of equations 
with added generalized quantum potential and account of a pressure term, for the oscillating wave packet. The 
quantum force applied on the fluid is calculated from a Gaussian fit of the density distribution. The figure gives 
the probability density in function of position (space grid from 90 to 140) and time (time steps from 1 to 32). 
In this simulation (near half a period), the amplitude of the pressure term is k p = 5. 



One finds that the addition of pressure leads to an oscillatory slight spreading of the wave 
packet, but that its main superfluid-like features are preserved, since it nearly recovers its 
shape after half a period. 




Fi gure 5: Result of the numerical integration of a Euler and continuity one-dimensional system of equations 
with added generalized quantum potential and account of a pressure term, for the oscillating wave packet. The 
quantum force applied on the fluid is calculated from a Gaussian fit of the density distribution. The figure gives 
the density distribution in function of the position (space grid form 120 to 160), for 12 equal time steps which 
cover a full period. In this simulation, the amplitude of the pressure term is k p = 1. One sees that the effect 
of pressure amounts to an oscillating spreading of the wave packet, which nearly recovers its shape after half a 
period. 



6 Full finite differences simulation 

The success of this first simple simulation leads us to attempt a more direct feedback in 
which the quantum force is computed by finite differences from the values of the density 
itself (while in the previous simulation we used an intermediate polynomial fit from which 
the force was analytically derived). 

To this purpose, we use a form of the generalized quantum potential and of the gen- 
eralized quantum force according to which they can be expressed in terms of only VlnP 



(or equivalently V In p) . Setting 

H = VlnP, 



(36) 



we find: 

Q = -V 2 (v.H + ^P 2 ) , (37) 

F Q = - VQ = V 2 [AH + (H.V)H] . (38) 

In one dimension it reads 

F *- v {-i^r + ^r^r)- (39> 

The numerical integration proceeds following the same lines as in the previous simulation, 
except for the first steps aiming at computing Fq , which are replaced by a finite difference 
calculation according to equation (38). Such a way to compute the force Fq to be applied 
on the fluid is therefore directly similar to its calculation in a real laboratory experiment 
from digitalized measurements of the density by a grid of detectors. Namely, we calculate 
successively, for all values of the position index j, 



2 5x 



H? = -^J"^- 1 , (40) 



then similar relations for positions Xj-i, Xj+i, Xj-2 and Xj + 2, then 

m-H 



then a similar relation for Q(xj+i, t n ), and finally 

The calculation of In p (from the continuity equation) and of V (from the Euler equa- 
tion) are the same as previously. We have attempted to use other more precise formulas 
for the calculation of the second and third order derivatives in the expression of Fq: this 
has led to essentially the same result. 

Despite, once again, the roughness of the chosen integration method, the result ob- 
tained is satisfactory, since the motion of a quantum non-spreading oscillating wave packet 
has been reproduced on about 1/4 of period before divergence due to the effect of com- 
puting errors (Figure 6). This result has been obtained without using the Schrodinger 
equation, but instead an apparently "classical" hydrodynamic Euler/continuity system 
with an externally applied generalized quantum potential. 

Adding a pressure term yields a similar result (i.e., reproduction of the motion of 
the wave packet on about 1/4 of period before divergence due to the effect of computing 
errors) which confirms the result obtained with the Gauss fitting method, namely, a partial 
oscillating spreading of the wave packet (Figure 7). 



Figure 6: Result of the numerical integration of a Euler + continuity one dimensional system with generalized 
quantum potential for the oscillating wave packet in an harmonic oscillator field. The quantum force applied 
on the fluid is directly calculated from the values of the density by finite differences. The density distribution 
obtained in function of position (space grid from 35 to 65) on about 1/4 of period has been used to reconstruct 
a full period (32 time steps). 
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Figure 7: Result of the numerical integration of a Euler + continuity one dimensional system with generalized 
quantum potential for the oscillating wave packet in an harmonic oscillator field and account of pressure (1/4 of 
period before stop due to computing errors). The quantum force applied on the fluid is here directly calculated 
from the values of the density by finite differences. The figure gives the probability density in function of space 
(grid from 100 to 130) and time (time steps from 1 to 13). A pressure term has been added (k p = 1, whose 
effect is a slight oscillatory spreading of the wave packet. 

7 Discussion and conclusion 

These preliminary simulations were intended to yield a first validation of the concept of a 
new kind of quantum-like macroscopic experiments based on the application to a classical 
system of a generalized quantum force through a retroaction loop [1]. They have given a 
positive results, since the expected quantum-type stable structure (here a non-spreading 
or slightly spreading oscillating wave packet) has been obtained during a reasonably long 
time of integration. These results, obtained by a rather rough integration method, are 
very encouraging since they give the hope that a real laboratory experiment should be 



possible to achieve. 

In the hydrodynamic case considered in this work, possible shortcomings are to be 
considered, such as the effects of finite compressibility, of vorticity, of viscosity at small 
scales, of density detector uncertainties, of the minimal time interval needed to perform 
the loop for the calculation and the application of the quantum force in a real experiment, 
etc... 

We have attempted here to have a first account of these uncertainties by taking a 
pressure term into account, by adding large random fluctuations in the initial conditions, 
then by adding large fluctuations at each time steps of the simulation. The results obtained 
were again encouraging, since, despite the pressure term and the large fluctuations, the 
overall coherence of the wave packet and its period were preserved. We shall in forthcoming 
works attempt to take into account these effects in more complete numerical simulations 
with improved integration schemes, apply the same general concept to other types of 
systems, then lead a real hydrodynamic laboratory experiment [12]. 

Provided such an actual experiment succeeds, it could lead to many new applications 
in several domains: didactic ones (teaching of quantum mechanics), laboratory physics 
(macroscopic models of quantum systems, simulations of atomic and molecular systems, 
study of the quantum to classical transition, laboratory astrophysics [13, 14], models 
of biological-like systems [15]), technology (development of a new devices having some 
quantum-like properties and behavior), self-organization (plama confinement, control of 
turbulence ?, etc.). 

Acknowledgements. The authors gratefully acknowledge very fruitful discussions with 
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